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A  MOLECULAR  DYNAMICS  INVESTIGATION  OF 
THE  STRUCTURAL  CHARACTERISTICS  OF 
AMORPHOUS  AND  ANNEALED  POLY (VINYLIDENE 
FLUORIDE)  AND  VINYLIDENE  FLUORIDE- 
TRIFLUOROETHYLENE  COPOLYMERS 


INTRODUCTION 

The  Navy  has  long  been  interested  in  the  polymers  poly(vinylidene  fluoride), 
PVDF,  and  poly (vinylidene  fluoride/ trifluoroethylene),  P(VDF-TrFE),  because  of 
their  piezoelectric  properties,  which  make  them  potential  candidate  materials  for 
active  and  passive  underwater  transducers.  Research  results  from  Naval  Undersea 
Warfare  Center  Division,  Detachment,  and  other  laboratories  have  demonstrated  that 
these  materials  hold  great  promise  with  their  prime  advantage  over  existing  piezo¬ 
electric  materials  being  their  low  mass,  high  processability,  and  low  cost.  However, 
experimentally  it  is  difficult  to  identify  and  optimize  the  crystallization  and  annealing 
conditions  utilized  in  the  construction  of  these  polymers. 

Within  this  laboratory  both  applied  and  theoretical  investigations  are  ongoing 
with  the  focus  towards  an  increased  understanding  of  these  materials.  The  goal  of 
the  project  reported  here  was  to  gain  insight  into  the  molecular  structure  of  the  two 
related  polymers.  Specific  interest  lays  in  establishing  structure/ property  relation¬ 
ships  to  ultimately  make  predictable  the  ideal  crystallization/ annealing  conditions  for 
specific  applications  as  a  function  of  polymer  composition. 

Details  of  polymer  processes  at  the  atomic  or  molecular  level  are  important  for 
a  more  thorough  understanding  of  bulk  material  properties.  However,  these  details 
are  not  directly  accessible  experimentally  because  the  specific  motions  of  atoms  in  a 
polymer  cannot  readily  be  evaluated  by  traditional  laboratory  techniques.  As  an 
alternative,  one  can  consider  computational  molecular  modeling.  The  theoretical 
method  of  choice  in  the  investigation  reported  upon  here  was  the  molecular  dynam¬ 
ics  (MD)  methodology.  Computer  simulations  employing  MD  allow  for  the  elucida¬ 
tion  of  dynamic  structural  change  of  polymers  at  the  molecular  level  as  a  function  of 
temperature.  In  this  communication,  we  report  on  the  MD  analysis  of  the  structures 
of  PVDF  and  P(VDF-TrFE)  as  a  function  of  temperature  and  relate  the  results  to  the 
crystallization  processes  observed  experimentally  for  each  of  the  polymers. 
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THEORY  AND  BACKGROUND 


NATURE  OF  THE  POLYMERS  INVESTIGATED 

The  polymer  PVDF  consists  of  many  repeating  units  of  the  monomer 
vinylidene  fluoride,  VDF.  When  VDF  monomer  units  are  mixed  and  reacted  with  a 
second  monomer  unit,  trifluoroethylene,  TrFE,  the  copolymer  P(VDF-TrFE)  is  formed. 
Numerous  copolymers  can  be  constructed  differing  from  one  another  in  the  relative 
proportions  of  the  two  monomers  as  well  as  their  connectivity.  The  composition  of 
the  copolymer  is  reported  as  the  percent  moles  (%mol)  of  vinylidene  fluoride.  For 
example,  75  mol%  P(VDF-TrFE)  designates  a  composition  with  a  VDF/ TrFE  mole 
ratio  of  3:1.  The  most  suitable  compositions  for  underwater  piezoelectric  applications 
are  those  with  VDF  levels  between  60  and  90  mol%,  as  has  been  shown  by  an  earlier 
molecular  modeling  investigation  in  this  laboratory.1,2  For  many  sonar  applications, 
P(VDF-TrFE)  copolymers  are  superior  to  PVDF  because  the  latter  requires  mechanical 
stretching  prior  to  poling,  adding  substantially  to  fabrication  costs.  Because  the 
copolymers  can  be  poled  without  stretching,  they  can  be  solvent-cast  or  injection- 
molded  into  desired  geometries  and  thicknesses  in  a  simple  and  low-cost  manner.3 


GENERAL  BACKGROUND  REGARDING  POLYMER 
STRUCTURE/PROPERTY  RELATIONSHIPS 

Long -chain  molecules  (polymers)  can  exist  in  two  general  states.  In  the  liquid 
state  there  is  normally  no  long  range  order,  and  chain  bonds  adopt  rotational  angles 
that  are  dictated  by  energetics  and  statistical  mechanics,  leading  to  the  evolution  of 
random  structures.4  The  centers  of  gravity  of  the  chains  themselves  are  also  random¬ 
ly  distributed  relative  to  one  another  .  However,  under  appropriate  conditions,  a 
chain  with  sufficient  periodicity  can  develop  order  along  a  portion  of  its  length.  An 
ensemble  of  such  chain  segments  can  then  order  itself  into  a  three  dimensional  array, 
and  the  system  is  said  to  have  crystallized.5  The  characteristic  macroscopic  conse¬ 
quence  of  a  bulk  polymer  possessing  considerable  liquid-like  or  amorphous  character 
is  rubber-like  elasticity.  In  contrast,  bulk  polymers  of  high  crystallinity  tend  not 
towards  elasticity  but  rigidity.  Many  organic  polymers  are  semi-crystalline,  pos¬ 
sessing  both  amorphous  and  crystalline  character,  and  macroscopic  properties  are 
intimately  related  to  the  degree  of  crystallinity. 

GENERAL  BACKGROUND  OF  POLYMER  STRUCTURAL  ANALYSIS 

The  spatial  arrangement  of  atoms  in  a  molecule,  designated  by  the  term 
configuration,  is  specified  in  terms  of  the  coordinates  within  three-dimensional  space 
for  each  atom.  The  coordinates  of  the  atoms  within  the  molecule  are  governed  by  the 


nature  of  the  atom-to-atom  connectivity,  or  chemical  bonding  between  atoms. 
Specifically,  the  lengths  of  bonds  and  the  angles  between  bonds  meeting  at  a  given 
atom  (bond  angles)  are  restricted  within  a  narrow  range.  The  potentials  affecting 
rotations  about  chemical  bonds  are  less  restrictive  than  those  governing  either  bond 
lengths  or  angles.  Large  molecules,  such  as  polymers,  owing  to  rotation  about 
bonds  of  the  molecular  framework,  possess  an  array  of  configurations,  or  conforma¬ 
tions.  The  term  conformation,  for  purposes  here,  relates  most  specifically  to  the 
torsional  (dihedral)  angles  of  the  polymer  chain.  The  polymer  can  be  described  as  a 
continual  array  of  dihedral  angles  going  from  one  end  of  a  chain  to  its  terminus.  (A 
dihedral  angle  is  the  bond  rotation  angle  around  the  axis  of  the  center  two  carbon 
atoms  between  two  planes,  each  defined  by  the  center  bond  and  one  terminal,  as  in 
C-C-C-C-)  One  conformation  for  a  given  polymer  may  have  repeating  C-C-C-C 
dihedral  angles  of  some  value  X,  while  a  second  conformation  for  the  same  polymer 
may  have  a  repeating  motif  with  a  dihedral  of  Y.  These  specific  conformations  of  the 
same  polymer,  called  polymorphs,  are  not  necessarily  of  equal  energy,  and  can 
interconvert  as  a  function  of  many  variables  (temperature,  elongation,  presence  or  ab¬ 
sence  of  solvent  or  plasticzers,  etc.).  Specific  polymorphs  for  a  given  chemical 
composition  are  typically  designated  with  Greek  prefixes  (i.e.,  the  PVDF  a-phase,  or 
PVDF  P-phase)  and  are  typically  further  specified  with  terminology  relating  to  the 
repeating  dihedrals.  The  specific  dihedral  nomenclature  pertinent  here  is  designated 
as  follows:  t  =  trans  or  -180°  dihedral,  and  g  -  gauche  or  ~45  -  90°  dihedrals. 

The  key  point  of  the  above  treatment  is  this:  from  the  viewpoint  of  funda¬ 
mental  polymer  chemistry,  the  macroscopic  properties  of  the  material  are  often  deter¬ 
mined  as  much  by  the  polymer  molecular  conformation  (i.e.,  polymorph)  as  its 
chemical  composition. 


DETAILED  STRUCTURAL  CHARACTERISTICS  OF  POLYMERS 
INVESTIGATED 

When  PVDF  is  either  crystallized  from  the  melt  or  solvent  cast,  the  predom¬ 
inant  polymorph  that  initially  forms  is  called  the  a-phase.6  The  conformation  of  the 
polymer  segments  within  the  crystalline  a-phase  unit  cell  is  tgtg',  with  the  chains  ori¬ 
ented  in  the  cell  such  that  no  net  polarization  exists  for  the  cell  (figure  1).  On  the 
basis  of  an  early  X-ray  crystallography  investigation,  Hasegawa  and  coworkers 
proposed  that  the  PVDF  a-phase  is  a  distorted  tgtg'  structure,  with  a  carbon  dihedral 
repeating  motif  of  179°  and  45°. 7  Doll,  Lando,  and  Bachmann  proposed  a  more 
"ideal"  tgtg'  conformation  with  g  and  g'  angles  of  ~60°.89  When  PVDF  is  uniaxially 
elongated,  the  a-phase  undergoes  polymorphism  to  the  P-phase.  Hasegawa  proposed 
that  the  P-phase  conformation  is  a  slightly  distorted  all- trans,  planar  zig-zag  with  an 
alternating  172°  and  -172°  repeating  motif.7  In  the  p-phase,  the  polymer  segments  are 
oriented  such  that  the  C-F  dipoles  are  additive,  resulting  in  a  polar  unit  cell  (figure 
1).  Hence,  to  become  piezoelectric,  PVDF  must  be  uniaxially  elongated  (to  induce  a- 
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to  (3-phase  polymorphism)  after  initial  crystallization  and  prior  to  electroprocessing. 

In  contrast,  the  copolymers  crystallize  directly  into  a  (3-phase  that  possesses  a  polar 
unit  cell.  The  copolymers  can  be  electroprocessed  into  a  piezoelectric  material 
immediately  after  crystallization. 

Several  theoretical  investigations  have  confirmed  the  existence  of  polymor¬ 
phism  in  PVDF.  Karasawa  and  Goddard,10  using  Hartree-Fock  calculations  and 
experimental  phonon  frequencies,  recently  reported  that  the  differences  in  energies 
between  four  polymorphs  of  PVDF  (including  the  a-  and  (3-phases)  are  within  1 
kcal/ mol/ monomer.  Farmer,  Hopfinger,  and  Lando,11  using  potential  energy 
calculations  to  determine  the  chain  conformation  and  packing  energies,  confirmed 
that  conversions  between  the  a-  and  (3-phases  should  be  energetically  feasible,  but 
that  increasing  the  concentration  of  head-to-head  defects  in  PVDF  should  favor 
formation  of  the  (3-phase. 

The  central  focus  of  our  ongoing  work  deals  with  the  structural  implications  of 
increasing  fluorine  content  in  P(VDF-TrFE)  copolymers.  The  fluorine  content  of  the 
copolymer  can  be  readily  controlled  by  altering  the  mol  percent  of  TrFE  with  respect 
to  VDF  in  the  generation  of  the  copolymer.  We  recently  reported  results  from 
molecular  mechanics  calculations  using  the  MM2  force  field  of  Allinger12  and  crystal 
packing  simulations  relating  to:  (1)  the  effect  of  fluorine  content  on  the  intramolecular 
energy  of  PVDF  and  various  copolymer  chain  segments;  and  (2)  the  relative  energies 
of  PVDF  and  copolymer  crystallites.2 

The  intramolecular  energies  for  polymer  chain  segments  were  found  to 
increase  with  increasing  fluorine  content  independent  of  conformation  (Flasegawa  or 
"ideal")  .  Interestingly,  for  PVDF,  the  Flasegawa  tgtg'  conformation  is  more  stable 
than  the  Hasegawa  and  "ideal"  dW-trans  structures.  In  contrast,  the  "ideal"  all -trans 
conformations  are  favored  for  each  of  the  copolymers.  The  crossover  point  corre¬ 
sponding  to  the  composition  where  the  Hasagawa  all-frans  conformation  becomes 
favored  over  the  tgtg'  conformation  is  ~90  mol%  VDF.  In  comparing  the  crystal 
packing  energies  for  the  Hasegawa  and  "ideal"  tgtg'  crystal  structures,  the  Hasegawa 
polymorphs  are  favored  for  all  compositions.  The  extent  to  which  the  Hasegawa 
tgtg'  is  preferred  decreases  somewhat  with  increasing  fluorine  content. 

The  salient  finding  from  our  previous  work  that  serves  as  the  motivation  for 
the  present  investigation  is  the  fact  that  the  PVDF  all- trans  conformation  (which,  as 
an  isolated  small  polymer,  is  less  stable  than  the  tgtg')  packs  more  favorably  into  a 
unit  cell  than  the  Hasegawa  tgtg'  conformation.  Thus,  it  may  be  that  the  all  -trans 
polymorph  is  the  global  energy  minimum  for  PVDF  even  though  the  polymer  is 
known  to  crystallize  from  the  melt  into  a  tgtg'  repeating  motif.  Thus,  it  is  reasonable 
to  assess  whether  intramolecular  interactions  exert  kinetic  control  on  the  crystal- 
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lization  behavior  of  the  polymer.  We  report  here  the  results  from  a  molecular 
dynamics  perspective  on  the  amorphous  and  the  annealed  structures  for  PVDF  and 
50  mol%  P(VDF-TrFE). 


EXPERIMENTAL 

THE  MOLECULAR  DYNAMICS  APPROACH 

In  a  molecular  dynamics  or  MD  simulation,  a  computer  algorithm  is  used  to 
calculate  the  future  positions  and  momenta  of  every  atom  in  a  molecule  based  on 
current  positions  and  momenta.13  The  positions  and  momenta  are  continually  propa¬ 
gated  in  a  series  of  equally-spaced  time  intervals  for  a  finite  time  chosen  by  the 
experimentalist.  In  this  fashion,  it  is  possible  to  simulate  conformational  change 
within  a  polymer  as  a  function  of  temperature  (i.e.  mimicking  crystallization  and/ or 
annealing)  in  a  dynamic  fashion  while  monitoring  potential  energy,  kinetic  energies, 
dihedrals,  chain  length  change  with  time,  etc. 

In  the  MD  approach,  Hamilton's  equations  of  motion  are  solved.  The  Simula-  " 
tion  first  determines  the  force  on  each  atom  (fj)  from  all  of  the  other  atoms  as  a 
function  of  time  at  a  certain  temperature: 


F  =-— 

'  *,* 


(i) 


where  V  is  a  potential  energy  function  and  rt  is  the  position  of  atom  i.  Then  the 
acceleration,  a,,  of  each  atom  is  determined  by  dividing  the  force  acting  on  it  by  the 
mass  of  the  atom,  mj( 


ar 


fi 


(2) 


The  approach  employed  in  the  present  work  is  called  the  leap-frog  algorithm. 
This  algorithm  updates  the  positions  of  atoms  and  the  velocities  for  the  next  step  by 
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the  following  equation: 


V;(f +-  A/) =v;(f-— Af) +Ata{t).  (^) 

2  2 

Let  v1/2  be  the  first  guess  of  a  velocity  at  t0-l/2At  and  x0  be  the  position  at  time  t0. 
The  "new"  velocity  v1/2  is  therefore 


Vl/2=V-l/2+^A^ 


(4) 


and  the  "new"  coordinate  x1  is 


*i=Vvi/2A'- 


(5) 


New  velocities  and  positions  can  be  specified  continuously  as  shown  below: 


V*2=VV2+aiA?> 

and 


(6) 


X2=Xl+V3/2 (?) 

At  each  integral  time,  the  new  coordinates  (represented  by  x  above)  and  the  accelera¬ 
tions  are  calculated;  at  half  integral  times,  the  velocities  are  determined.  The 
molecular  dynamics  trajectory  can  be  visualized  as  a  set  of  atoms  with  specific 
coordinates  and  velocities  evolving  according  to  equations  (3)  to  (7). 

The  kinetic  energy  (K)  is  then  defined  in  terms  of  the  velocities  and  masses  of 
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the  atoms,  i.e.. 


K=-mvf.  (*) 

2  ‘ 

The  kinetic  energy  at  time  t  is  calculated  using  the  current  velocity,  which  is  approxi¬ 
mated  by 


(9) 


The  total  energy  of  the  system,  or  Hamiltonian,  is  the  sum  of  the  kinetic  and  poten¬ 
tial  energies: 


H(rj>)=K(p)+V(p),  (1U) 

where  r  is  a  set  of  Cartesian  coordinates  and  p  is  the  momenta  of  the  atoms. 

In  a  typical  molecular  dynamics  simulation,  the  positions  and  velocities  are 
assigned  to  the  atoms.  The  velocities  reflect  the  target  temperature  for  the  simula¬ 
tion.  Equilibration  proceeds  with  partition  between  potential  and  kinetic  energy  and 
a  thermal  distribution  of  velocities.  Equilibration  may  involve  thousands  of  time 
steps.  In  the  production  part  of  the  run,  as  time  evolves,  the  coordinates  and 
velocities  at  each  time  step  are  stored.  This  part  of  the  analysis  allows  the  study  of 
the  properties  and  geometries  of  new  structures  as  they  emerge  with  time. 


COMPUTATIONAL  METHODOLOGY 

All  computations  were  performed  on  a  Silicon  Graphics  4D/ 35  IRIS 
workstation  utilizing  commercial  software  produced  by  Hypercube  (HYPERCHEM) 
and  Cerius.  The  MD  computations  were  performed  at  the  MM+  level  utilizing  the 
force  field  of  Burkert  and  Allinger.14  All  structures  were  evaluated  atomistically. 


SIMULATION  CONDITIONS 

Two  criteria  were  utilized  in  determining  the  optimum  temperature  for  the 
simulations  of  the  amorphous  polymers.  First,  the  dihedral  angles  of  the  structure, 
upon  heating,  must  not  exhibit  "memory  effects"  with  respect  to  the  input  structure 
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(i.e.,  all  of  the  dihedral  angles  within  the  structure  must  change  markedly  from  the 
input  dihedrals  as  evidenced  by  real-time  observation  of  the  dynamic  structure  as  it 
is  heated).  Second,  heating  the  structure  slightly  above  or  slightly  below  the 
chosen  temperature  should  not  change  the  distribution  of  dihedral  angles  significant¬ 
ly  from  those  attained  at  the  chosen  temperature.  Both  of  these  criteria  were  chosen 
in  order  to  simulate,  as  closely  as  possible,  the  melt  conditions  of  the  polymer.  The 
heating  temperatures  utilized  do  not  quantitatively  relate  to  their  real-world  analogs 
(i.e.,  an  input  temperature  of  1800  K  does  not  relate  at  all  to  a  real-world  tempera¬ 
ture  of  1800  K),  but  rather  represent  values  utilized  by  the  program  to  manipulate  the 
extent  of  atom  displacement  per  cycle.  At  a  heating  temperature  of  2200  K,  the  chain 
structure  begins  to  distort  (i.e.,  begins  to  change  bond  lengths  and/ or  angles  in  an 
unrealistic  fashion),  while  at  temperatures  below  1650  K  memory  effects  seem  to  be 
present.  At  temperatures  between  1700  and  2100  K,  the  trends  in  the  distribution  of 
dihedral  angles  are  similar,  although  not  identical.  The  relative  trends  in  structure  as 
a  function  of  temperature  (from  hot  to  cold),  whether  from  2000  K  to  270  K,  or  from 
1750  K  to  270  K,  do  not  qualitatively  change.  In  view  of  these  considerations,  the 
temperature  chosen  for  all  simulations  was  2000  K.  The  1  picosecond  (ps)  heating 
time  was  chosen  for  expediency  as  longer  heating  times  produced  similar  results, 
since  longer  times  use  considerable  computer  time.  Once  the  polymer  was  heated,  it 
was  allowed  to  equilibrate  for  3  ps.  The  3  ps  equilibration  period  was  chosen 
because,  based  on  these  observations,  it  was  well  into  the  region  in  which  the  kinetic, 
potential,  and  total  energies  of  the  system  were  constant  (1  ps  would  have  been  suffi¬ 
cient  based  upon  this  criteron).  It  should  be  noted  that  the  choice  of  a  trajectory 
length  of  the  order  of  picoseconds  imposes  a  limitation  on  a  molecular  dynamics 
simulation,  in  that  only  molecular  motions  occurring  on  a  short  time-scale  4,  such  as 
bond  stretching  vibrations,  are  sampled.  Only  changes  in  atomic  positions  are 
monitored  in  the  picosecond  regime. 

In  simulating  the  annealing  of  the  polymers,  a  10  ps  cooling  time  was  utilized 
in  going  from  2000  K  to  270  K.  This  represents  a  rather  slow  annealing  and  is 
typical  of  simulations  of  this  nature. 


RESULTS  AND  DISCUSSION 

MOLECULAR  DYNAMICS  OF  RELATED  DIMERS 

As  a  preamble  to  evaluating  polymers,  a  simplified  MD  investigation  of  the 
dimers  of  vinylidene  fluoride  and  vinylidene  fluoride/  trifluoroethylene  was  under¬ 
taken.  Each  dimer  initially  possessed  a  180  degree  C-C-C-C  dihedral  angle.  The 
input  structures  were  heated  to  high  temperature  (2000  K)  over  a  1  ps  time  frame 
and  allowed  to  equilibrate  at  that  temperature  for  an  additional  3  ps,  at  which  point 
the  dihedral  angle  was  evaluated.  Then  the  structure  was  allowed  to  anneal  for  a  10 
ps  time  frame  to  300  K,  and  again  the  dihedral  angle  was  evaluated.  From  table  1,  it 
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is  clear  that  the  hot  vinylidene  fluoride  dimer  possesses  a  gauche  conformation  which 
upon  annealing  remains  gauche,  whereas  the  hot  vinylidene  fluoride/ trifluoroethylene 
dimer  possesses  a  gauche  conformation  which  upon  cooling  becomes  trans. 


MOLECULAR  DYNAMICS  OF  SINGLE-CHAIN  POLYMERS 

In  an  effort  to  model  a  small  PVDF  polymer  whose  structure  is  reflective  of 
the  amorphous  polymer  (i.e.,  from  the  melt),  a  128  methylene-unit  long  end-capped 
polymer  was  generated  (initially  with  an  all-trans  structure)  and  heated  to  high 
temperature  (2000  K)  over  a  1  ps  time  frame  and  allowed  to  equilibrate  at  that 
temperature  for  an  additional  3  ps.  In  the  evaluation  of  the  resulting  structure,  the 
C-C-C-C  dihedral  angles  were  analyzed  (as  absolute  values  such  that  both  -45°  and 
45°  are  reported  as  45°,  170°  and  -170°  reported  as  170°,  etc.).  The  results  (table  2) 
reveal  a  preference  for  trans- like  dihedrals  (between  145°  and  180°),  although  the 
structure  is  also  rich  in  gauche-like  bonds  (dihedrals  between  37°  and  90°).  The 
gauche/trans,  or  g/t  ratio  for  the  amorphous  polymer  is  0.73,  suggesting  that  the 
polymer  at  high  temperature  has  substantial  gauche-character. 

When  the  above  polymer  is  allowed  to  anneal  over  a  long  time  frame  (10  ps),  „ 
it  adopts  a  new  conformational  structure.  The  analysis  of  the  dihedrals  indicates 
that  the  polymer  becomes  more  ordered  upon  cooling  in  that  the  percent  of  dihedrals 
other  than  those  that  are  either  trans-  or  gauche- like  is  reduced  (from  24%  to  2%). 

The  percent  of  contribution  from  both  trans-  and  gauche-like  bonds  increases  with 
cooling,  though  the  magnitude  of  the  increase  for  gauche-like  bonds  is  greater  (table 
2).  Thus,  the  g/t  ratio  of  the  annealed  structure  is  increased  upon  cooling  ( g/t  =  0.82 
for  the  annealed  and  0.73  for  the  hot  polymer).  In  figure  2,  a  histogram  comparing 
the  dihedrals  associated  with  both  the  hot  and  cold  PVDF  and  the  trends  leading  to 
greater  order  and  enrichment  in  gauche-character  with  cooling  are  graphically  repre¬ 
sented. 

An  identical  approach  was  used  in  the  analysis  of  50  mol%  P(VDF-TrFE),  the 
results  of  which  are  shown  in  figure  3  and  table  2.  Once  again,  the  polymer 
becomes  much  more  ordered  upon  cooling.  Flowever,  unlike  PVDF,  the  percent  trans 
character  increases  dramatically  upon  cooling  but  the  percent  gauche  character 
remains  constant  (table  2).  Thus,  the  resulting  g/t  ratios  indicate  that  the  copolymer 
is  enriched  in  trans- character  upon  cooling  (g/t  for  amorphous  structure  is  0.64  while 
that  for  the  annealed  structure  is  0.46). 

In  comparing  hot  PVDF  to  hot  P(VDF-TrFE),  the  PVDF  is  only  slightly  richer 
in  gauche  conformations  (figure  4).  In  figure  5,  the  comparison  between  the  two 
polymers  after  annealing  reveals  the  same  trend  (that  PVDF  is  richer  in  gauche- 
character).  On  cooling,  the  homopolymer  becomes  richer  in  gauche-character,  whereas 
cooling  the  copolymer  does  not  produce  any  additional  gauche  bonds.  These  results. 
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which  are  in  accord  with  the  experimental  trends  observed  upon  the  cooling  of  both 
PVDF  and  50  mol%  P(VDF-TrFE),6  are  probably  due  to  the  differences  in  bond,  angle, 
torsional,  van  der  Waal,  and  electrostatic  energies  between  PVDF  and  the  more 
fluorinated  copolymers. 


MOLECULAR  DYNAMICS  OF  MODIFIED 
SINGLE-CHAIN  POLYMERS 

A  modified  starting  structure  for  PVDF,  one  with  a  tgtg'  repeating  motif,  was 
utilized  to  compare  with  the  results  from  the  earlier  all-frans  PVDF.  This  structure 
was  chosen  for  investigation  because  the  tgtg'  polymorph  for  PVDF  is  well  known.7 
When  the  tgtg'  PVDF  structure  is  first  heated  and  then  annealed  under  the  identical 
conditions  as  those  stated  earlier,  the  gjt  ratio  increases  from  0.69  to  0.90  (amorphous 
and  annealed,  respectively,  in  agreement  with  results  from  the  earlier  experiment 
utilizing  the  all-trans  PVDF  input  structure). 

A  second  type  of  modified  input  structure  was  also  utilized,  based  upon  an 
insight  gained  on  having  evaluated  the  aforementioned  single-chain  polymers.  In 
each  of  the  single-chain  structures  investigated,  the  polymers  bend  back  upon 
themselves  a  number  of  times  such  that  small  regions  of  the  chain  are  interacting 
with  one  another.  The  extent  of  chain  folding  is  small  and  difficult  to  quantitate  for 
128  methylene-unit  polymers.  However,  were  it  feasible  to  evaluate  much  longer 
polymers,  it  is  likely  that  the  effect  would  be  more  pronounced.  Because  investigat¬ 
ing  larger  polymers  is  impractical  (because  computational  times  measured  in  weeks 
would  be  required),  an  alternative  input  structure  was  utilized.  The  structure 
possesses  two  parallel  and  dynamic  chain  segments  (~7  A  apart)  where  one  end  of 
each  chain  segment  is  tethered  to  a  common,  fixed  U-shaped  link.  This  structure 
makes  possible  a  simulation  that  includes  enhanced  intermolecular  interaction 
because  the  two  chain  segments  are  parallel,  in  close  proximity  to  one  another,  and 
dynamic. 

The  analysis  of  the  U-shaped  PVDF  polymer  reveals  that  the  structure  becomes 
more  ordered  upon  cooling  with  an  increase  in  both  the  percent  contribution  from 
trans-  and  gauche-like  bonds  (table  3).  The  g/t  ratio  for  amorphous  "folded  over"  U- 
shaped  PVDF  is  0.64,  which ,  upon  annealing,  remains  relatively  constant  (changing 
only  slightly  to  0.65).  The  analysis  of  the  U-shaped  50  mol%  copolymer  again 
reveals  increased  order  upon  cooling,  with  the  g/t  ratio  remaining  constant  at  0.53 
(the  structures  themselves  change,  though  the  statistical  ratio  remains  the  same). 

That  PVDF  possesses  greater  gauche-character  in  both  the  amorphous  and  annealed 
structures  than  does  the  copolymer  is  consistent  with  both  the  single  chain  computa¬ 
tions  and  the  experimental  behavior  of  the  true  polymers.6 
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ELONGATION  OF  tgtg'  PVDF 


When  PVDF  is  cast  as  a  thin  film,  it  initially  possesses  a  tgtg'  repeating  motif, 
consistent  with  both  the  results  of  our  earlier  investigation2  regarding  crystal  packing 
and  the  MD  results  of  this  study.  When  the  thin  film  is  mechanically  stretched,  the 
initial  tgtg'  structure  undergoes  polymorphism  to  an  all -trans  structure.6  To  model 
such  a  phenomenon,  a  tgtg'  PVDF  tetramer  was  constructed  with  initial  dihedral 
angles  at  45°,  180°,  -45°,  180°  (corresponding  to  the  Hasagawa  distorted  tgtg'  geome¬ 
try  experimentally  observed  for  PVDF  crystallites7  based  upon  x-ray  defraction  and 
consistent  with  our  earlier  crystal  packing  computations2).  All  bond  angles  and 
distances  were  then  allowed  to  relax  via  geometry  optimization  (with  fixed 
dihedrals).  The  procedure  for  geometry  optimization  is  discussed  in  an  earlier 
report.1  The  resulting  structure  possessed  an  end-to-end  distance  of  9.252  A.  This 
initial  structure  was  then  incrementally  elongated  in  0.25  A  steps  to  an  ultimate  end- 
to-end  distance  of  10.252  A.  With  each  elongation,  the  complete  structure  was 
allowed  to  relax  (with  the  end-to-end  distance  fixed),  mimicking  the  elongation  of  the 
crystallite.  Throughout  the  elongation,  all  the  C-C-C-C  dihedral  angles  were 
monitored.  From  table  4,  a  clear  and  systematic  polymorphism  from  the  tgtg' 
structure  to  an  all-trans-like  structure  is  evident  as  each  gauche  bond  smoothly 
converts  to  a  trans  bond  upon  elongation. 


CONCLUSIONS 

In  an  earlier  study,  we  used  molecular  modeling  calculations  to  show  that 
energetics  plays  a  major  role  in  determining  the  conformational  nature  of  the 
predominant  polymorph  as  PVDF  or  its  copolymer  with  trifluoroethylene  is  annealed 
from  the  melt  or  solution.1 2  Although  calculations  of  intramolecular  steric  energies 
were  consistent  with  the  crystallization  behavior  of  PVDF  and  its  copolymer,  crystal 
packing  energies  suggested  that  the  all  -trans  polymorph  is  the  global  minimum  struc¬ 
ture  for  PVDF  and  P(VDF-TrFE).  The  lack  of  agreement  between  the  experimental 
observation  that  PVDF  is  favored  to  crystallize  directly  into  the  a-phase  rather  than 
the  {3-phase  suggested  to  us  the  possibility  of  a  kinetic  barrier  hindering  the  forma¬ 
tion  of  the  more  thermodynamically-stable  (3-phase.  The  MD  simulations  reported 
here  confirm  that  crystallization  and  crystal  packing  of  PVDF  and  the  copolymers  is 
kinetically  controlled.  Theoretically,  these  calculations  suggest  that  the  PVDF  a-phase 
is  indeed  a  metastable  state,  and  that  the  barriers  for  formation  of  the  all  -trans  are 
much  higher  for  PVDF  than  for  the  copolymers  (figure  6).  The  kinetic  barriers  for 
formation  of.  this  state  probably  consist  of  unfavorable  torsional  conformations.  On 
the  practical  side,  this  result  explains  why  mechanical  stretching  along  the  polymer 
backbone  is  a  requirement  for  formation  of  the  piezoelectric  (3-phase  of  PVDF  and  not 
for  the  copolymers.  As  summarized  in  figure  6,  the  pathway  to  the  PVDF  (3-phase 
must  proceed  directly  through  the  a-phase,  whereas  for  the  copolymers  direct 
crystallization  into  the  [3-phase  is  the  favored  route.  The  nature  of  the  conformational 
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barriers  separating  the  a-  and  (3-phases  in  PVDF  has  yet  to  be  elucidated.  However, 
as  suggested  by  Flory  in  his  classic  work,4  very  small  differences  in  the  energies  of 
gauche  and  trans  conformations  in  fluoroalkanes  may  lead  to  rather  dramatic  differ¬ 
ences  in  their  physical  properties.  No  doubt  these  small  differences  play  a  vital  role 
in  the  annealing  and  crystallization  behavior  of  PVDF  and  P(VDF-TrFE). 
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DRAWING 


Figure  1.  Structures  ofPVDF  a-  and  /3-Phase  Polymorphs 


#  of  Dihedrals 


Hot  PVDF  vs. Cool  PVDF 


Figure  2.  Histogram  Showing  Distribution  of  Dihedrals  of  Hot  and  Cold  PVDF 
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#  of  Dihedrals 


Hot  50%  vs. Cool  50% 


Dihedral  Angle 

Figure  3.  Histogram  Showing  Distribution  of  Dihedrals  of  Hot  and 
Cold  50  mol%  P(VDF-TrFE)  Copolymer 
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#  of  Dihedrals 


Hot  PVDF  vs.  Hot  50% 


Dihedral  Angle 

Figure  4.  Histogram  Showing  Distribution  of  Dihedrals  of  Hot  PVDF 
vs.  Hot  50  mol%  P(VDF-TrFE) 
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it  of  Dihedrals 


Cool  PVDF  vs.  Cool  50% 


Dihedral  Angle 


Figure  5.  Histogram  Showing  Distribution  of  Dihedrals  After  Annealing 
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Table  1.  Calculated  Dihedral  Angles  of  Simple  Dimers 


Dimer 

Dihedral  angle,  deg 

Vinylidene  fluoride  (hot) 

-67.83 

Vinylidene  fluoride  (cool) 

-50.96 

Vinylidene  fluoride/trifluoroethylene  (hot) 

68.29 

Vinylidene  fluoride/trifluoroethylene  (cool) 

171.53 

Table  2.  Calculated  gauche-  and  trans-Dihedral  Angles  of  128-Methylene  Unit 

Single  Chain  Polymers 


Polymer 

%t 

%g 

%  other 

all- trails  PVDF  (hot) 

0.73 

44 

32 

24 

all- trails  PVDF  (cold) 

0.82 

54 

44 

2 

all- trails  50  mol%  P(VDF- 
TrFE)  (hot) 

0.64 

47 

30 

23 

all -trans  50  mol%  P(VDF- 
TrFE)  (cold) 

0.46 

67 

30 

3 

tgtg'  PVDF  (hot) 

0.69 

43 

30 

27 

tgtg'  PVDF  (cold) 

0.90 

51 

47 

2 
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Table  3.  Calculated  gauche-  and  trans-Dihedral  Angles  of  U-Shaped  Polymers 


Polymer 

8/* 

%t 

%g 

%  other 

all -trans  PVDF  (hot) 

0.64 

47 

30 

23 

all-trans  PVDF  (cold) 

0.65 

58 

38 

4 

all-fnms  50  mol%  P(VDF- 
TrFE)  (hot) 

0.53 

51 

27 

22 

all -trans  50  mol%  P(VDF- 
TrFE)  (cold) 

0.53 

63 

34 

3 

Table  4.  Elongation  of  a  tgtg' PVDF  Tetramer  as  Monitored  by 
Geometry  Optimization 


End-to-end 

o 

distance  (A) 

C-C-C-C  dihedral  angle,  deg 

1 

2 

3 

4 

5 

9.252 

45 

180 

-45 

180 

45 

9.502 

83 

175 

-84 

-174 

83 

9.752 

98 

173 

-91 

-172 

89 

10.002 

154 

171 

-97 

-168 

124 

10.252 

170 

169 

-168 

-166 

176 
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